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ABSTRACT 

We investigate the global nonaxisymmetric Rossby vortex instability in a 
differentially rotating, compressible magnetized accretion disk with radial density 
structures. Equilibrium magnetic fields are assumed to have only the toroidal 
component. Using linear theory analysis, we show that the density structure 
can be unstable to nonaxisymmetric modes. We find that, for the magnetic field 
profiles we have studied, magnetic fields always provide a stabilizing effect to the 
unstable Rossby vortex instability modes. We discuss the physical mechanism of 
this stabilizing effect. The threshold and properties of the unstable modes are 
also discussed in detail. In addition, we present linear stability results for the 
global magnetorotational instability when the disk is compressible. 

Subject headings: accretion, accretion disks — magnetohydrodynamics — insta- 
bilities — waves 



1. INTRODUCTION 

Radial density structures can be common in astrophysical disks. For example, it is 
thought that a transition region exists in protostellar disks that can be caused by different 
degrees of coupling of the magnetic field with the disk material (Gammie 1996; Terquem 
2008). Tidal interactions between protoplanets and protoplanetary disks can also give rise 
to dips/gaps in protoplanetary disks (Goldreich & Tremaine 1980; Lin & Papaloizou 1986; 
Ward 1997). At the inner edge of black hole accretion disk, the sharp density contrast 
between the plunge region and the accretion disk can also be treated as a density transition 
structure. 
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Rossby wave instability or Rossby vortex instability (RVI) in thin Keplerian disks with 
density structures in the hydrodynamic limit has been studied in the linear theory (Lovelace 
et al. 1999; Li et al. 2000). The existence of unstable modes has been found to be associated 
with the radial gradients of an entropy-modified potential vorticity profile. RVI in disks will 
form vortices and shocks in the nonlinear limit (Li et al. 2001). Recently, RVI associated 
with a dip/gap is studied by de Val-Borro et al. (2007). The intrinsic mechanism of the 
instability is the corotation amplification caused by over- reflection (Goldreich et al. 1986). 
Unlike Papaloizou & Pringle (1984) instability, the RVI does not depend on the reflection 
boundary conditions (see Li et al. 2000 for a detailed discussion). 

Magnetic fields are supposed to be present in accretion disks and would greatly change 
the dynamics of disks. Previous magnetized disk instability analyses, both local and global, 
have mainly considered smooth disks without structures (e.g., Balbus & Hawley 1991; Ogilvie 
& Pringle 1996). Global magnetorotational instabihty (MRI) with vertical and azimuthal 
magnetic fields have been studied by Curry & Pudritz (1996) and Ogilvie & Pringle (1996), 
respectively. They both found that global MRI are localized and Curry & Pudritz (1996) 
gave an explanation in terms of the Alfven resonance positions in the disk. Curry & Pudritz 
(1996) also found that global MRI is very sensitive to the boundary conditions. Recently Pino 
& Mahajan (2009) propose to reduce the dependence on boundary condtions by restricting 
the rotation rate change only in a narrow range in radius. Magnetized disk instability with 
step-shaped density structure has been considered by Tagger & Pellat (1999), but they just 
considered purely poloidal magnetic fields. This kind of accretion-ejection instability only 
occurs in strongly magnetized disks with plasma /3 < 1. The growth rates are typically 
small, ~ several percent of the Keplerian frequency. Due to the disk differential rotation, 
the dominant component of magnetic fields should be toroidal (Balbus & Hawley 1998). How 
would the toroidal magnetic fields change the behavior of the RVI is still an open question. 

In this paper, we make such an attempt to study the magnetic field effects on structured 
disks and extend the RVI study into disks with toroidal magnetic fields. The paper is 
organized as follows. In §2 we describe the equilibrium disk with density structures. In §3, we 
give the linear analysis of the magnetized RVI. We also discuss the effects of compressibility 
on global MRI. Conclusions are given in §4. 

2. Disk Equilibrium and Structure Profile 

We use a cylindrical (r, 0, z) coordinate system. The equilibrium disk is axisymmetric 
and in steady state, with unperturbed velocity vq = v^e^. The vertical stratification is 
neglected and is assumed to be uniform. Only the azimuthal component of the magnetic 
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fields is present. For the axisymmetric equilibrium disk, the radial force balance reads 
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where $ is the gravitational potential of the central object and the disk self-gravity is not 
considered in the present paper. Here p, p and are mass density, gas pressure and 
toroidal magnetic field, respectively. We define VLq to be the Keplerian angular velocity at 
ro and Vq = ro^o. The length, time and mass density scale are normalized by tq, fig ^ and 
Po, respectively. The disk is assumed to be isothermal with a constant temperature. The 
isothermal sound speed is chosen as Cg, and Pq = pocl- The toroidal magnetic field is taken 
as 

= ASor-'^ , (2) 

where Bq = y/po- The strength and gradient of the toroidal magnetic field is controlled by 
A and a, respectively. We typically take a = 1 in our paper. At the characteristic radius ro, 
the plasma P in our formulation is approximately P ~ Pq/{\^BI/2) = 2/A^, but, in general, 
/3 is a function of radius. 

We focus in this paper on the configurations of disks with density dip/gap with a toroidal 
magnetic field and study the stability of such equihbrium disks. For simplicity, we model 
the density dip/gap with a Gaussian profile 



p = Po < 1 - - 1) exp 



(3) 



where G (1 < G < 2) and A specify the amplitude and width of the density structure, re- 
spectively. The above profile is motivated by the fact that the protoplanet's tidal interaction 
with the gaseous disk would induce a gap in the disk. The perturbation by a protoplanet 
leads to the excitation of spiral density waves at Lindblad resonances, which carry an an- 
gular momentum flux. The waves deposit the angular momentum flux when the waves are 
dissipated. As a result the outer disk receives angular mometum from the protoplanet and 
the inner disk loses angular momentum to the protoplanet. The outer (inner) disk gas gains 
(loses) angular momentum and moves outward (inward). When the tidal torque is greater 
than the viscous torque of the disk, a surface density dip/gap would be formed in the vicinity 
of the protoplanet. In numerical simulations the gap profile is similar to a Gaussian profile. 
Such configurations could also be relevant to the radial border of active zone and dead zone, 
where strong coupling of the magnetic field and disk material gives rise to bigger accretion 
rate in the active zone compared to the dead zone. As a result, there would be a density 
jump at the boundary between the active zone and the dead zone. The instability behavior 
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of the density jump at this border is quite similar to the outer edge of the density gap. For 
simphcity, we just consider the gap situation in this paper. Note that the instabihty behavior 
of the outer edge of the density dip/gap is very simihar to the density jump case considered 
in Li et al. (2000). Wc take G = 1.5, Tq = 1, A/tq = 0.05 and Cg = 0.07. The inner and outer 
radii of the disk arc taken as Tin — 0.4 and Tout — 1.6, respectively. In Figure 1, we give the 
equilibrium disk with a dip/gap. The four panels are P/Po{ro), fl{r)/flK{r), ij-) / VL\;{r-) 
and /Tj^I, respectively, where k is the radial epicyclic frequency and — The 
quantity «;^/(En) is the profile of potential vorticity (PV). We can observe that there are 
two minima at r ~ 0.9 and 1.1 in the PV profile related to the inner and outer edges of the 
gap, which imply that there may exist two unstable regions associated with both the inner 
and outer edges. 



3. Linear Analysis of Vortex Instability 

Small perturbations to the inviscid compressible Euler equation are considered. The 
mass density is p — Pq -\- 5p, the gas pressure is p = po + ^Pi the velocity is v = Vq + 
u — Vlre^f, + u and B = B^e^f, + 5B, where 5p, 5p, u and 5B are perturbed mass density, 
gas pressure, velocity and magnetic field, respectively. The magnetohydrodynamic (MHD) 
equations of the compressible disk are 

^ + pV-v = 0, (4) 

= -Ivp- V$ + i(V X B) X B , (5) 
Dt p p 

— = V X (v X B) . (6) 

Note that the factor Air is absorbed in the definition of magnetic field. 

We linearize the above equations by taking perturbations proportional to exp(im0 + 
ikzZ — icut), where m — ±1, ±2, ... and are the azimuthal and axial wavenumber, respec- 
tively, and cij = a;^ + ^7 is the mode eigenfrequency. We use ^ = {Sp + B^SB^) /p as the 
perturbed total pressure divided by density {SB^ is the perturbed azimuthal magnetic field) 
and as the Lagrangian displacement. The linearized perturbation equations for and ^ 
read 

^ = Au^r + Au^ , (7) 

= A2lCr + >l22* . (8) 
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Fig. 1. — The equilibrium structure of the disk with a dip/gap. In this example the solid 
lines correspond to A = and a = 0, the dashed line correspond to A = 2.0 and a = 2. 
Usually, we choose a = 1, but when a = 1, the dashed lines will overlap with solid lines. So 
we choose a = 2 to make the differences appreciable. From top to bottom, the panels are 
gas pressure, P/Pq, angular velocity, fl/flx, square of epicyclic frequency /Vt\-, potential 
vorticity K^/iT.iX). In the first panel, the solid line and dashed line overlap. 



The four coefficients in the above two equations are 
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^a,2 — ^^C^ 



Note that in the axisymmetric incompressible disks with uniform density background, the 
quantity ^l^ j. reduces to — -^(S^/r)^, which is the same as given by Chandrasekhar (1961). 
When taking the incompressible limit, equations (JTj) and ([H]) reduce to (2.13) and (2.14) in 
Ogilvie & Pringle (1996). The two equations (JTj) and ([H]) can be combined to get a single 
second order differential equation with respect to \E', 



^ + B{r)^ +C(r)^ = , (15) 

where 

D/'«,\ I A I /, 
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Cir) = - (auA2i - AnA22 + A'^2 - 

\ ^21 

and the prime denotes the derivative with respect to r. 



3.1. Axisymmetric Stability 

Before we proceed to investigate the vortex instability produced by inflexion points in 
the PV profile, we need to ensure the equilibrium is stable to axisymmetric perturbations. 
To this end, we use the sufficient condition for local stability, i.e., the generalized Rayleigh 
criterion 

K\r) + nl^;,o. (16) 

This criterion can be readily seen in the incompressible axisymmetric uniform density back- 
ground model. In such a case the two coupled equations become (we temporarily keep in 
this subsection in order to see the generalized Rayleigh criterion and suppress kz in the later 
RVI analysis) 

D,^r = and = {a^ ^1,.% , (17) 

where = D + ^ and D = Eliminating \E' from the above two equations, we arrive at 

{DD^-kl% = -'^,{n' + nlr)^r (18) 

Simple variational principle analysis gives that a necessary and sufficient condition for ax- 
isymmetric perturbations to be stable is that + be positive throughout the whole 
radius range (e.g., Chandrasekhar 1961). In all of our following examples, we require that 
the above asxisymmetric stablility criterion is met. We will see that the stronger the mag- 
netic field, the bigger the quantity Q"^ ^ and it is this term that stabilizes the hydrodynamic 
RVI. 
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3.2. Method of Solving Linear Eigenvalue Problem 

Since the eigenfrequency uj is in general complex, equations ([7]) and ([8]) are a pair 
of first order differential equation with complex coefficients which are functions of r. We 
choose to use the relaxation method to solve these equations (Press et al. 1992). In this 
method, the ODEs are replaced by finite difference equations on a mesh of points covering 
the domain of interest. The relaxation method needs an initial trial solution that can be 
improved by a Newton-Raphson technique. Iterations are carried out by carefully designed 
Gaussian elimination adapted to block diagonal matrix. After iterations the initial trial 
solution will gradually converge to the two point boundary eigenvalue problem. For the 
MRI calculation, we use rigid boundary conditions. For the RVI calculation, the boundary 
conditions implemented are such that waves propagate away from the density structure in 
both inner and outer parts of the disk (e.g., Li et al. 2000). 

3.3. 3D Results of Linear Analysis on MRI 

We solved equations ([7]) and ^ in two cases. One corresponds to a three dimensional 
infinite cylinder with fc^ 7^ and the other corresponds to two dimensional thin disk with 
kz = 0. In this section we focus on the MRI and we treat the disk density as uniform (i.e. 
G = 1) and kz 7^ 0. Only the global nonaxisymmetric MRI is considered. The difference 
between our model and Ogilvie & Pringle (1996) (hereafter OP) is that they considered the 
incompressbile limit while our model is a compressible one. 

The equilibrium setup given by OP has, in normalized units, po = = r^/^, = 
A^r^^, As is a constant and taken as 0.2. Here As is related to A in equation (2) by As = A Cg. 
We take k^ = 28 and m = 5 in the following analysis for we can get relatively higher growth 
rate with this choice of parameters. In order to get to the realistic value of sound speed 
for disks, such as Cg = 0.07, we first recover OP's results by taking Cg = 100 and A^ = 0.2. 
Then we gradually reduce the sound speed from Cg = 100 to Cg = 0.5 and keep Xb = 0.2. 
We find that when the sound speed reaches Cg = 0.5, we can not find unstable modes any 
more. From a physical point of view, this means that the magnetic fields become too strong 
and the MRI is suppressed. Then when we keep Cg = 0.5 and gradually reduce A^ from 0.2 
to 0.1, we can find unstable MRI again. After this, we keep Xb = 0.1 and gradually reduce 
the sound speed Cg from 0.5 to 0.2. The MRI unstable modes growth rates decrease as we 
reduce the sound speed Cg from 0.5 to 0.2. Repeating the above process for several times, 
we gradually get to a realistic value of sound speed Cg = 0.07. In Figure [21 we show how 
the growth rate of MRI varies with the sound speed while the magnetic field strengh is fixed 
with Xb = 0.08. 
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Once we find the unstable mode with the realistic sound speed for disks, we can slowly 
vary the magentic field strength to find the dependence of growth rates on magnetic field. 
Figure [3] gives the effects of the magnetic field strength on the growth rate for a particular 
mode with = 28 and m = 5 for Cg = 0.07. Most unstable modes are achieved with 
moderate magnetic field strength. This is different from the behavior of RVI, whose growth 
rate is a monotonically decreasing function of the magnetic field strength (see §4.3.2). And 
the peaks are achieved at different magnetic field strength for the inner and outer modes. 
The inner mode reaches the greatest growth rate around A = 0.4 or = 0.028 and the 
outer mode around A = 1.71 or Xb = 0.12. Surprisingly, we find that the plasma (3 < 1 for 
the high A situation in Figure [31 Usually MRI is suppressed in low jS magnetized gas, here 
we can see that gas compressibility extends the MRI to the low (3 regime. 

The eigenfunctions of unstable modes are quite similar to the results of Ogilvie & Pringle 
(1996), the compressible global MRI modes are also localized at the boundaries, which has 
been explained by Curry & Pudritz (1996) in terms of Alfven resonance positions in the 
disk. Figure H] shows examples of the MRI eigenfunction of the above mentioned two most 
ustable modes. We find that the amplitudes concentrate either at the inner or at the outer 
boundary (that is why these modes are named). The most unstable inner mode frequency 
is ujr/i^^o) = 3.8077, the radial postions where u — mfl{r) = ±k^Ca{r) are ri = 0.3979 
and r2 = 0.4219 and the corotaion radius is at Vd = 0.4102. We note that ri is outside the 
computaional domain, and the amplitude is mainly confined between the inner boundary 
Tin = 0.4 and r2. The most unstable outer modes frequency is uJr/ (^^o) = 0.5124. The radial 
positions where uj — mf2(r) = ±A;0Ca(r) are = 1.4564 and r4 = 1.6573 and the corotaion 
radius is at = 1.5617. We note that, in this case, r4 is outside the computational domain. 
The amplitude is mainly confined between and the outer boundary rout = 1-6. These 
observations are essentially the same as the results of Curry & Pudritz (1996). We expect 
that the highly localized MRI would not affect the RVI much even though it has greater 
growth rates than RVI, because the radial density structure is far from the disk boundaries. 
We have also investigated the azimuthal wave number dependence of the MRI growth rate. 
Figure [5] gives the variation of growth rate with the azimuthal wave number m of these two 
modes. The inner m = 3 mode and outer m = 6 mode have the greatest growth rates. We 
use the rigid boundary conditions = at both boundaries. We also tried the outflow 
boundary conditions used in the RVI calculation, but we can not find unstable modes any 
more. So global MRI is quite sensitive to boundary conditions as shown by other studies as 
well (e.g.. Curry & Pudritz 1996). 
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2. — Dependence of growth rate of MRI on the compressbihty for a fixed magnetic field 
- 0.08. kz — 28 and m — 5 are chosen. 




3. — Dependence of the MRI growth rate 7 on the magnetic field strength for the 
5 modes with = 28.0, Cg = 0.07. 
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Fig. 4. — The left (right) two panels show the inner (outer) mode eigenfunction and \E'. 
The real part is shown by the solid line, the imaginary part by dashed line. The inner mode 
is mainly between inner boundary r^^ = 0.4 and r2 = 0.42 while the outer mode is between 
ri = 1.46 and the outer boundary rout = 1-6. 
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Fig. 5. — Dependence of the MRI growth rate on the azimuthal wave number m. Solid line 
is for the inner mode, dashed line for the outer mode. 
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3.4. 2D Results of Linear Analysis on Rossby Vortex Instability 

For disks with radial density structures, we restrict tlie analysis in the 2D (r, 0) plane 
{kz = 0). In this limit, MRl is not excited. For the following analysis, we will focus on RVI 
with magnetic fields. 

3.4-1- Representative hydrodynamic examples 

We first show the hydrodynamic results for a representative case with G = 1.5, Cg = 0.07 
and A/ro = 0.05. Solving equations ([7]) and ([8]), we find many unstable modes. Figure [6] 
shows the dependence of mode frequency Ur and growth rate 7 on the azimuthal mode 
number m for the inner edge mode and the outer edge mode, respectively. The overall 
behaviour of the hydrodynamic instability is in good agreement with results of Li et al. 
(2000) and de Val-Borro et al. (2007). 

3-4-2- Results of RVI with Magnetic Fields 

With the inclusion of magnetic fields, we calculate the growth rate of the unstable 
modes. In Figure [7] and El we present the results of the effects of different field strengths on 
the frequency and growth rate for different azimuthal wavenumber modes associated with the 
inner and outer edge of the dip/gap, respectively. The variation of the frequency and growth 
rate of the m = 5 outer edge mode with different magnetic field strength is presented in Figure 
[9l For this particular mode, both the frequency and growth rate decrease monotonically with 
the increase of magnetic field. When the magnetic field is strong enough, we can see that 
the Rossby vortex instability will be completely suppressed by the presence of magnetic 
field. From Figure [9l when the parameter A ~ 0.88, or when the plasma (3 is approximately 
2/A^ ~ 2.6, the RVI is almost completely suppressed. 

The eigenfunctions of the m = 5 unstable modes with A = 0.1 are shown in Figure [TU] 
and Figure [11] for the inner edge mode and outer edge mode, respectively. The unstable 
inner edge mode has a growth rate ~ 0.1181 and a real frequency uJr/{mVLQ) k, 1.1287. 
The unstable outer edge mode has a growth rate 7/f2o ~ 0.1211 and a real frequency 
UrlirnVto) ^ 0.8730. We have used outward-propagating sound wave boundary conditions 
to obtain these eigenfunctions. The relative phase shift between real and imaginary parts 
indicates this propagation. Two dimensional distribution of the perturbation are shown in 
Figure [121 We can identify that vortices develop around the radial density structure. The 
azimuthal pressure gradient is crucial for the formation of anticyclonic vortices (see Fig. 5 in 
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Li et al. (2001) for a detailed explanation). When the equilibrium azimuthal magnetic fields 
are present, magnetic fields restrict the development of radial velocity thus the formation of 
the vortices. As a result the growth rate is reduced by the azimuthal magnetic fields. 

3.4-3. Stabilizing effect of Magnetic Fields and Physical Mechanism for the Instability 

As we mention before in the axisymmetric analysis, the quantity ^ plays an essential 
role in stabilizing the RVI. Note that Q"^ ^ can be written as: 

Usually, the instability takes place near the density structure in the disk, where the density 
scale length is much shorter than the disk radius. So in the above equation, the dominant 
terms are those related to the density length scale. We can see that 

""..-*5e^^|p|j(^)(^)>0. (20, 

From this eqution we can see that, when the instability ocurrs at either the density decreasing 
inner edge or the density increasing outer edge, the contribution of Q^ j. is always positive, 
that is, it will decrease the instability growth rate and stabilize the RVI. 

In Figure [121 we give both and ^ as a function of radius, we can see that stronger 
magnetic fields give larger value of fl'^ ^ which stabilizes the RVI. To investigate the magnetic 
field gradient on the behavior of the instability, we also tried different values of a = —1,0,1,2, 
the stabilizing effect of toroidal magnetic fields still holds. 

The basic physical mechanism for these modes is essentially the same as the hydrody- 
namic RVI (Lovelace 1999; Li et al. 2000). Although in the nonaxisymmetric analysis, the 
"potential" — C(r) becomes complex, it is still of guidance to plot the real part of the "po- 
tential", C(r) in equation f|T5l) . The real part of the function — C(r) for different magnetic 
field strengths is shown in Figure O For the outer edge mode, the negative "potential well" 
around r/ro = 1.05 is the unstable region. When an unstable mode is excited in the "poten- 
tial well" of the unstable region, the two positive potential peaks at two sides of the unstable 
region cause this mode to be evanescent in such regions. These two trapping regions would 
partially act as refiection boundary for the amplification to work, although there will be a 
finite probability for the mode to tunnel through the "potential barriers" (as the "potential" 
becomes again negative when it goes far from the unstable region). 

The "potential well" of —Cir) in the wave equation can also confirm stabilizing effects 
of magnetic fields. Comparison of Figure [13] and Figure [TH shows that the increase of ^ 
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at r = 1.05 causes the potential to become narrower and shallower. This means that the 
excitation of the instability is weaker. Interestingly, another potential well seems to have 
developed at r = 1.1, which may mean that there exists another branch of modes in our 
equations. But we suspect that this branch of modes is stable and therefore is not the focus 
of this paper. 

3.4-4- Instability Threshold 

Consider now the dependence of the growth rate on the dip/gap amplitude G- As the 
amplitude G decreases, the growth rate of the instability is expected to decrease. Figure [13 
shows the calculation of the growth rate and mode frequency of hydrodynamic Rossby vortex 
mode as a function of G for the m = 5 unstable mode. The threshold value is Gth = 1-32 
and 1.26 for the inner edge and the outer edge mode, respectively. The threshold of the 
outer edge mode is slightly lower than that of the inner edge mode. Note that the value of 
Gth depends on A and c^, which are not discussed here. For A = 0.4, Gth = 1-38 and 1.28 for 
the inner edge and the outer edge mode, respectively. When magnetic fields are included, 
Gth increases. 

4. Conclusions 

We have carried out a linear analysis of the magnetized RVI associated with an axisym- 
metric, local radial density structure in a thin accretion disk. The flow is made unstable due 
to the existence of local extreme in the radial profile of the potential vorticity. Depending on 
the parameters, the unstable modes are found to have substantial growth rates ~ 0.1f2(ro), 
where tq is the location of the surface density dip/gap. We also find that strong toroidal 
magnetic fields would stabilize the RVI induced by the radial density structures in disk. In 
situations where only weak magnetic fields exist, such as protoplanetry disks, the existence 
of the magnetized RVI and their nonlinear outcome as vortices (e.g., Li et al. 2005) indicate 
that the flows in the co-orbital region are more complicated than what is expected from the 
linear analysis. If the disk plasma f3 is usually large (a few tens), the corresponding A is 
about 0.3, from Figure [HI we can see that the instability is only slightly depressed by the 
toroidal magnetic fields and the RVI is still active even with the toroidal magnetic fields. 
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Fig. 6. — Dependence of the RVI mode frequency cUr and growth rate 7 on the azimuthal 
mode number m. The sohd and dashed hne are for the inner edge (A) and outer edge mode 
(B), respectively. is zero in this case. 
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Fig. 7. — Dependence of the mode frequency uOr and growth rate 7 on the azimuthal mode 
number m for the inner edge mode with different strength of magnetic field, a is set to be 
1. The grow rate decreases with the strength of magnetic field, but the frequency increases 
slightly with the strength of the magnetic field. For small magnetic field A = 0.01, it is 
indistinguishable from the nonmagnetic case. 
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Fig. 8. — Same as Fig. [7] but for the outer edge mode. 
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Fig. 9. — Dependence of the m — 5 outer edge mode frequency ujr and growth rate 7 on 
the strength of magnetic field. As the magnetic field increases, the growth rate diminishes. 
When A ~ 0.88, or plasma /3 ~ ~ 2.6, the RVI is suppressed. 
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Fig. 10. — Eigenfunction for the inner edge mode (Mode A) of RVI with magnetic fields 
A = 0.1. Shown are the perturbed density, the radial and azimuthal magnetic perturbations 
and the perturbation function ^ for ni = 5. The dashed line is the real part, the dot-dashed 
line is the imaginary part, and the solid line is the amplitude. 
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Fig. 11. — Same as Fig. [TO] but for the outer edge mode. 
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Fig. 12. — Two dimensional distribution of the gas pressure perturbation based on linear 
theory of the outer edge m = 5 mode with A = 0.1 and a = 1.0. Arrows indicate the 
perturbation magnetic vector near tq. Units are arbitrary. 
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Fig. 13. — The radial profile of /t^ and Q'^^ of the m = 5 outer edge mode for different 
magnetic field strength. As the magnetic field increases, ^2^^, increases around where the 
excitation of instabihty takes place (both the inner edge and outer edge) and stabilizes the 
RVI. 
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Fig. 14. — The "potential well" around r = 1.05 becomes narrower and shallower with the 
increase of magnetic field strength. The excitation of the RVI becomes weaker. For small 
magnetic field A = 0.01, it is indistinguishable from the nonmagnetic case. 
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Fig. 15. — Dependence of the mode growth rates on the amphtude of the surface density 
gap/dip G for the m = 5 unstable modes. The vanishing of the growth rate for G < Gth 
indicates the threshold for the RVI. When magnetic fields are included, the threshold vlaues 
of G increase. 
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